g'9/2 nuclei and neutron-proton interaction 
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We have performed shell-model calculations for nuclei below ^'"^Sn, focusing attention on the two 
N = Z nuclei ®^Cd and ^^Pd, the latter having been recently the subject of great experimental 
and theoretical interest. We have considered nuclei for which the Ogg/2 orbit plays a dominant role 
and employed a realistic low-momentum two-body effective interaction derived from the CD-Bonn 
nucleon-nucleon potential. This implies that no phenomenological input enters our effective Hamil- 
tonian. The calculated results for ^^Pd are in very good agreement with the available experimental 
data, which gives confidence in our predictions for Cd. An analysis of the wave functions of both 
^^Cd and ^^Pd is performed to investigate the role of the isoscalar spin-aligned coupling. 



PACS numbers: 21.60.Cs, 21.30.Fe, 27.60.-t-j 



I. INTRODUCTION 

In a recent work [l| three excited states in the N — Z 
nucleus ^^Pd, lying quite far from the stability line, were 
observed. This is a remarkable achievement since the 
N — Z nuclei play a special role in understanding the nu- 
clear effective interaction, in particular the interplay and 
competition between the isovector and isoscalar compo- 
nents. In fact, in this case neutrons and protons occupy 
the same orbitals, which gives rise to a large overlap of 
their wave functions. In this context, the main question, 
which is still a matter of debate [3j Q ' i^ ^^'^ existence 
of strongly correlated T — np pairs, similarly to the 
well-known case of neutron and proton pairs coupled to 
J = and T = 1. 

In Rcf. 1] it was pointed out that the main fea- 
ture of the measured levels of ^^Pd is their approximate 
equidistance. An interpretation was given of this feature 
in terms of a shell-model calculation with an empirical 
Hamiltonian. This reproduces very well the experimental 
levels and yields wave functions built manly from J = 9 
np pairs, which has been seen as evidence for a spin- 
aligned np paired phase. 

The results of 0| have immediately attracted much 
attention and the role of the isoscalar np pairs in the 
low-energy structure of the N = Z nuclei close to doubly 
magic ^'^'^Sn has been investigated in some very recent 



theoretical papers |4-|6|. These studies confirm substan- 
tially the dominance of J = 9, T = pairs in the low- 
lying yrast states of the N — Z nuclei with four, six and 
eight holes below ^""^Sn. 

The relevance of the isoscalar component of the n 



np 



interaction has been stressed in a subsequent work 
where the 16"*" "spin-gap" isomer in ®^Cd has been identi- 
fied and its origin explained as due to the strong influence 
of this component. 

Another interesting outcome of the above works is that 
in all nuclei above *^Ru the low-lying yrast states can be 
essentially described in terms of the single Ogg/2 shell. 
This situation is of course reminiscent of the so-called 



/7/2 nuclei that have been the subject of a large num- 
ber of theoretical studies ever since the seminal paper by 
McCuUen, Bayman and Zamick f8|. It seems therefore 
appropriate to call these special nuclei below ^'^'^Sn "39/2 
nuclei" . 

Some ten years ago, we performed J9| shell-model cal- 
culations for A^ = 50 nuclei immediately below ^"^Sn 
employing a realistic effective interaction derived from 
the Bonn-A nucleon-nucleon (NN) potential by means 
of a G-matrix formalism. In that work we took as model 
space for the valence proton holes the four levels f^/2, 
P3/2, Pi/2 and 39/2 of the 28-50 shell. Our results turned 
out to be in very good agreement with the available ex- 
perimental data. Since then, however, a new paradigm 
for realistic shell-model calculations has been developed 
which consists in renormalizing the strong short-range 
repulsion of the bare NN potential through the so-called 
Mow-k approach ^10|J. Furthermore, high-precision po- 
tentials have been constructed which fit the pp and np 
scattering data with a x^ /datum w 1. 

The exciting new findings mentioned before have stim- 
ulated us to perform modern realistic shell-model calcu- 
lations for nuclei below ^"^Sn, with particular attention 
focused on ^^Pd and on the heavier N = Z nucleus ^®Cd, 
for which theoretical predictions are likely to be verified 
in the not too distant future. Based on the dominant 
role of the 59/2 orbit in the low-lying states of the nuclei 
considered in the present study, we have restricted our 
calculations to the single 59/2 shell. This permits a more 
transparent analysis of the structure of the wave func- 
tions, especially for ^^Cd, in terms of either the [nn](8'[pp] 
or [np] (g) [np] coupling schemes. The comparison between 
these two approaches is instrumental to understand the 
role of the J = 9 np pairs. It may also be worth recalling 
that within the fpg space the choice of the single-hole 
energies is not an easy task @, since there is no spec- 
troscopic information on the single-hole nuclei ^^In and 
99Sn. 

In Sec. II we focus attention on the gg/2 effective in- 
teraction employed in our calculations. Results for ^^Cd 



and ^^Pd are presented in Sec. Ill, where we also give a 
detailed analysis of the structure of wave functions. Sec. 
IV provides a summary and concluding remarks. 



II. 



59/2 



EFFECTIVE INTERACTION 



We assume ^°"Sn as a closed core and let the neutron 
and proton holes move in the single gg/2 orbit. Our two- 
body effective interaction is derived within the frame- 
work of the time-dependent degenerate linked-diagram 
perturbation theory [il[ starting from the high-precision 
CD-Bonn NN potential [12|. This potential, that as all 
modern NN potentials contains high-momentum non- 
perturbative modes, is renormalized by constructing a 
low-momentum potential Viow-k- This is achieved by 
integrating out the high momentum modes of the bare 
potential down to a cutoff momentum A = 2.1 fm~^. 
Then the smooth V^ow-k potential plus the Coulomb force 
for protons is used to calculate the two-body matrix ele- 
ments of the effective interaction by means of the Q-box 
folded-diagram expansion [11| , with the Q-box including 
all diagrams up to third order in Viow-k- These dia- 
grams are computed using the harmonic oscillator basis 
and considering intermediate states composed of all pos- 
sible hole states and particle states restricted to the six 
proton and neutron shells above the Fermi surface. The 
oscillator parameter hw is 8.55 MeV, as obtained from 
the expression huj = AbA~'^/^ - 25A'^/^ for A = 100. 
The shell-model calculations have been performed by the 
NUSHELLX code [H. 



TABLE I: Proton-proton, neutron-neutron, and proton- 
neutron matrix elements of 14ff in the 39/2 orbit (in MeV). 
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-0.271 
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0.317 


0.083 


0.066 
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-0.404 
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0.459 


0.221 


0.210 
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-1.402 



We show in Table H] the two-body matrix elements of 
the effective interaction. Owing to the Coulomb force, 
there is no isospin symmetry, the nn matrix elements be- 
ing more attractive than the pp ones by about 250-350 
keV, which agrees quite well with the results of previous 
works [IJ, [15 where the effective interaction was deter- 
mined by a least-squares fit to experimental energies. As 
regards the T = 1 np matrix elements, they differ only by 
at most 100 keV from the nn ones. This reflects the fact 
that for np and nn valence holes/particles the Coulomb 
force acts only through Q-box diagrams starting at sec- 



ond and third order, respectively. From Table U we also 
see that the most attractive matrix elements correspond 
to the J = 0, J = 1, and J — 9 states, the J = matrix 
element being the largest one while the other two have 
about the same magnitude. This dominance of the J = 
matrix element stems from the fact that our interaction 
is derived for the single 179/2 orbit. The same feature is 
shown by the three different 179/2 interactions of Ref. [j| 
while this is not the case for interactions defined in the 
fpg space, as for instance that developed in |16| . 




FIG. 1: Calculated proton hole-neutron hole multiplet in ^*In. 



The np matrix elements, relative to the J = energy, 
correspond to the excitation spectrum of ^*In, which we 
find interesting to show in Fig. [TJ This has unfortunately 
no experimental counterpart, so cannot be used to test 
our interaction. We see, however, that the np multiplet 
exhibits a downward parabolic behavior, which is typical 
of nuclei with one proton-one neutron valence particles 
or holes, as for instance '*^Sc or ^^Co with two particles 
and holes, respectively, in the /7/2 orbit. The observed 
multiplet in these nuclei shows the same behavior as that 
we have found for ^*In, the only difference being a larger 
dispersion in the energy values which is due to the greater 
attractiveness of the np interaction in lighter systems. 




FIG. 2: (Color on line) Experimental and calculated proton 
hole-neutron particle multiplet in ^''Nb. 



For a test of our np matrix elements, we compare in 
Fig. [5] the observed multiplet in ^"Nb [I3| with the cal- 
culated one. The latter is, in fact, directly related to the 



matrix elements of the hole-hole np interaction through 
the Pandya transformation [1^. Note that several 1+ 
states have been observed in ^"Nb at low excitation en- 
ergy. Following the suggestion of Ref. [4)1 , we report the 
fifth one at 2.126 MeV while exclude the observed 0+ 
state at 5.008 MeV, which is too high in energy to make 
trustable our interpretation in terms of the single gg/2 
model. Fig. [5] shows that the agreement between the- 
ory and experiment is very good, the calculated energies 
overestimating the experimental values by at most 150 
keV in the case of the 2+ state. The latter, however, 
is likely to be admixed with configurations outside our 
model space. 




FIG. 3: Experimental and calculated spectra of the A'' = 50 
isotones '^^Cd, ^'^Ag, and ^«Pd. 



We now focus attention on the pp matrix elements by 
comparing the spectra of the three TV = 50 isotones ^®Cd, 
^^Ag, and ^^Pd with the experimental ones. We include 
all observed levels for ^^Cd while positive-parity yrast 
levels below 4.5 MeV for the two latter nuclei. This is 
done in Fig. [31 where we see that the agreement between 
theory and experiment is quite satisfactory up to about 
3.7 MeV. We overestimate the energy of the first excited 
state in the three nuclei by less than 100 keV while the 
energies of all other levels, except the 12+ in ^^Pd which 
is predicted at more than 400 keV below the experimental 
one, are underestimated by an amount which increases 
when the mass number decreases reaching at most 250 
keV. In concluding this section, it is worth noting that 
the accuracy of the present gg/2 calculation is similar to 
that of [i|, where the spectrum of ^^Pd up to the 10+ 
state was calculated in the fpg model space. 



III. RESULTS 

The results obtained for ^°Nb and iV = 50 isotones 
have given us confidence in our effective interaction, at 
least as regards its predictive power for low energy states. 
We have then performed calculations for the two N — Z 
nuclei ^^Cd and ^^Pd. The calculated spectra are re- 
ported in Figs, m and [5^, respectively, together with a 
comparison with the experimental data of [l|| for ^^Pd. 
In view of the recent experimental finding [7|] we have 



reported yrast states in ^^Cd up to I^ = 16+, while for 
^^Pd we have included two more states with respect to 
those identified in [l|. 
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FIG. 4: Calculated spectrum of ''"^Cd. 
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FIG. 5: Experimental spectrum of ^^Pd compared with the 
results of calculations with: a) all matrix elements, see Ta- 
ble U b) nn, pp, np T — 1 matrix elements; c) nn, pp, np 
T = matrix elements. 



The three observed states in ^ Pd are very well re- 



produced by the theory. Of course, one may not expect 
the same kind of agreement for high-energy states. It is 
worth mentioning, however, tliat our results account for 
the isomeric nature of the 16+ state identified in ^^Cd. 
More precisely, we find that the location of this state be- 
low the I^ ~ 12+ and 14+ states can be traced to the 
isoscalar np component of the interaction, in agreement 
with the results of [Ij. In particular, we have verified 
that a strong attractive J = 9 matrix element is essential 
to make the 16+ state isomeric. This is related to the 
fact that for the 12+, 14+, and 16+ states the average 
number of isoscalar J — 9 pairs (see Sec. III.B for the 
definition of this quantity) increases with angular mo- 
mentum, reaching its maximum value for the 16+ state. 
This is, in fact, 2.5 to be compared with 2.2 and 2.0 for 
the 14+ and 12+ states, respectively. Therefore a suffi- 
ciently attractive J — 9 matrix element pushes the 16+ 
state below the two lower-spin states. Relevant to this 
discussion is the fact that the structure of the wave func- 
tions for the 12+, 14+, and 16+ states, and consequently 
the average number of isoscalar J = 9 pairs, does not 
depend on the size of the J = 9 matrix element. 

Let us now focus on states up to I^ = 10+ for ^^Cd 
and ^^Pd. For both nuclei, we find an almost regularly 
spaced level sequence up to /'^ = 6+, with a slightly 
reduced 6+ — 4+ spacing. Then, the 8+ — 6+ spacing be- 
comes even smaller while the 10+ — 8+ one increases con- 
siderably. These features, which up to the 6+ state find 
a correspondence in the experimental data for ^^Pd, are 
only slightly more pronounced in ^^Cd. This similarity 
may be seen as an indication that the same correlations 
come into play in their low-energy spectra. We shall dis- 
cuss in more detail the results for both nuclei in the two 
next subsections. 



"Cd 



We start with ^^Cd, for which a simple and clear anal- 
ysis of the wave functions can be performed using the 
orthogonal basis formed by products of nn and pp states. 
The structure of the wave functions in terms of this basis 
set is shown in Table |lll while their overlaps with the 
[{np)9{np)9]i state are given in Table link . 



TABLE II: Overlap of the calculated /'' = 0+, 2+, 4+, 6+, 8+, 
and 10+ yrast states in ^^Cd with the [{nn)Jn{pp)Jp]i states, 
expressed in percentage. Only components with percentage 
> 10 are reported. 

('■'■ri, Jp) 
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TABLE III: Overlap of the calculated J'' = 0+, 2+, 4+, 6+, 
8+, and 10+ yrast states in Cd with the [(np)9(np)9]/ state, 
expressed in percentage, obtained using: a) Vij{np) of Tab. [H 
b) one-half the original value of V9(np); c) twice the original 
value of V9(np). 
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Note that for each angular momentum I an orthonor- 
mal basis formed by products of two np-pair vectors can 
been constructed, the [{np)9{np)9\i state being one of 
them. These vectors are simply related to the [nn] ® [pp] 
basis through 
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i I jI )\{nn)J,,{pp)Jp-I) , (1) 

^n ^p -* 

where [J] — (2 J + 1) and N denotes the normalization 
factor. From Table |TT] we see that the wave functions 
of ^^Cd expressed in terms of the [nn] ® [pp] basis are 
strongly fragmented. In fact, large seniority-4 compo- 
nents are present in the ground as well as in the excited 
/'^ — 2+, 4+, 6+, and 8+ states. More precisely, their 
percentage in the ground state is 43% while in the other 
states is not less than 34%. The 10+ state is character- 
ized by an admixture of different seniority-4 components, 
each with a weight not exceeding 18%. On the other 
hand. Table llllh shows that, when written in terms of 
an [np] ® [np\ basis, the ground and the first two ex- 
cited states are largely dominated by the [{np)9{np)9\i 
component. The weight of this component, however, is 
significantly smaller for the other three states having the 
minimum value (6%) for the 8+ state. These results for 
^^Cd are in line with those of Ref. [4i]. 

The large fragmentation evidenced in Table [Tl] is of 
course related to the np interaction and in this connection 
it is interesting to find out what is the role of the J — 9 
matrix element VQ{np). To this end, we have redone our 
calculations with two different values of VQ{np)^ namely 
increasing and reducing VQ{np) by a factor of 2. A similar 
analysis was done in 0, 0] ■ 

We find that for the reduced value of VQj{np) the 6+ 
and 8+ states come down in energy getting close to the 
4+ state. As a matter of fact, in this case the group 
of the 4+, 6+, and 8+ levels concentrates in a small en- 
ergy range separated by a large energy gap from both 
the 2+ and 10+ states. The spectrum of ^^Cd (up to 
/■^ = 8+) becomes then similar to that of ^*Cd, the wave 



functions being, however, still significantly seniority ad- 
mixed. The weight of the seniority-4 components is no 
smaller than 29%. On the contrary, a doubled value of 
VQ{np) leads to almost equidistant levels, with an energy 
separation of about 1 MeV, the only exception being the 
8+ — 6+ spacing which is ~ 300 keV smaller. This evo- 
lution toward an equidistant-level spectrum comes along 
with a larger fragmentation of the wave functions. The 
only state which has a less admixed nature is the 10+ 
state. We find that the percentage of the [(nn)0(pp)0]o 
component in the ground state reduces to 51%, while 
that of the seniority-two components, [{nn)I{pp)0]i and 
[{nn)0{pp)I]i, in the I = 2+, 4+, 6+, and 8+ states 
ranges from 33% to 63%. These values, however, re- 
main significantly large showing that pairing is still in 
the game. 

We now examine the infiuence of Vginp) on the dom- 
inance of isoscalar J = 9 pairs in the wave functions of 
^^Cd. To this end, we have also reported in Table IIIII 
the overlaps with the [{np)9{np)9]i state obtained us- 
ing one- half (b) and twice (c) the original Vg{np) value. 
We see that a larger value of Vg (np) leads to an increase 
of the overlap for all the states. However the overlap 
for the 8+ state does not go beyond 27%. Moreover, it 
should be noted that for the three lowest-lying states the 
overlaps do not become significantly smaller even when 
Vg (np) is reduced. This is related to the structure of the 
[{np)9{np)9]i^o,2,4 states in terms of the [nn] (g) [pp] basis 
[see Eq. (1)]. Therefore, as already pointed out in [J|, 
both dynamics and geometry are crucial to the presence 
of J = 9 pairs in the wave functions of ^^Cd. 



^Pd 



The simple analysis done for ^^Cd cannot be performed 
for ^^Pd with four neutron and four proton holes. In this 
case, we shall first discuss the effects of the T — and 
T — 1 components of the interaction on the calculated 
spectrum. This study was already performed by Ceder- 
wall et al. [1|, and we have found it worth verifying if a 
realistic effective interaction confirms their results. 

In Fig. [5] the results of the full-interaction calculation 
(a) are compared with those obtained by removing sepa- 
rately the T = (b) and T = 1 (c) np matrix elements. 
We see that in case (b) the excited states up to I'^ — 8+ 
are compressed in a smaller energy interval, about 1 MeV 
to be compared with 2 MeV of the full calculation. Ac- 
tually, the spectrum of ^^Pd evolves toward that of the 
neutron closed shell nucleus, as it was the case for ^^Cd 
when using a reduced value of Vg{np). On the other 
hand, when we exclude the T — 1 np matrix elements all 
the excited levels move down, but the spectrum keeps the 
same structure as that obtained from the full calculation. 

Our findings are in line with those of Ref. [l|], con- 
firming the more relevant role of the T = versus the 
T = 1 np component. More specifically, we have verified 
that the addition of the sole J = 9 np matrix element to 



the interaction between identical particles is sufficient to 
produce a spectrum very similar to that of Fig. [St. We 
may therefore conclude that the structural makeup of the 
^^Pd spectrum is mainly determined by the combined ac- 
tion of J = 9 np, nn and pp matrix elements. Needless 
to say, a quite distorted, highly compressed spectrum 
would result from ignoring the interaction between iden- 
tical particles. 

In this context we have tried to better understand how 
different isoscalar and isovector pairs contribute to pro- 
duce the spectrum of Fig. [5^. To this end, we have used 
the relation 



^i = J2^C'j{np)Vjinp)+C'jipp)Vj{pp)+C'jinn)Vjinn)] , 
J 

(2) 
where the energy of a given state is written in terms 
of the average numbers of nn, pp, and np pairs, 
C5(ij)'s, defined as C'j{ij) =< ^/ (^^Pd) [ [(a^ a])j x 
{aiaj)j]o I V'/ (^^Pd) >. In Eq. (2) the matrix elements 
in the three different channels appear explicitly, since, as 
mentioned in Sec. II, our effective interaction includes 
the Coulomb force. For the sake of simplicity, in the fol- 
lowing we shall not distinguish between nn, pp, and np 
isovector pairs and take as Vj, with J even, the mean 
value of the three corresponding matrix elements. 




FIG. 6: (Color on line) Average number of the isoscalar (a) 
and isovector (b) {gQ/2)^J pairs, Cj, as a function of the an- 
gular momentum I of the yrast states in Pd. 



In Figs. [6^ and |6}d, we show the average number of 
the isoscalar and isovector pairs for the six considered 
yrast states. We see that the curves corresponding to 
J = 8 and 9 lie significantly higher than the others, as 
was already observed in [5|. We draw attention here on 
the fact that the curves of Fig. [6^ are almost flat when 
compared to those of Fig. [B)d. This means [see Eq. (2)] 
that the isovector pairs contribute to the level spacings 
more significantly than the isoscalar ones. In particular, 
a main role is played by the J = and 2 pairs, whose 
corresponding matrix elements are much larger than the 
other ones. 



TABLE IV; Contributions (in MeV) of isoscalar and isovec- 
tor pairs to the energy level spacings of ®^Pd. Values < 0.13 
are omitted. Column S gives the sum of contributions while 
column AE the energy spacing obtained from the full calcu- 
lations. See text for details. 
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0.13 
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0.66 


8+ -6+ 


-0.64 
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0.18 


0.56 


0.52 


0.47 


10+ - 8+ 


1.42 


-0.26 


0.13 




-0.42 


0.88 


1.00 



The contributions to the five spacings arising from the 
different pairs are reported in Table IIVI vifhere we only 
include values larger than 130 keV. The two last columns 
allows to compare the calculated spacings of Fig. [S^ with 
those obtained by summing the contributions reported 
in the table. We see that they do not differ significantly. 
From Table lTVl it appears that the most important contri- 
butions to the energy spacings arise from the J = pairs, 
although those from J = 2 and 9 cannot be ignored at 
all, being particularly relevant for the 6 — 4 and 8 — 6 
spacings. In this regard, it should be kept in mind that 
the structure of the wave functions, and consequently the 
specific action of the interaction in the J = and 2 chan- 
nels, is strongly infiuenced by the size of the J = 9 matrix 
element. 

In concluding this discussion, we should mention that 
the content of isoscalar spin-aligned pairs in the wave 
functions of ^^Pd is still significant even when Vg{np) 
is suppressed. This, as already emphasized for ^^Cd, is 
related to geometrical features. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this work, we have performed a shell-model study 
oi N — Z nuclei below ^°°Sn that can be described in 
terms of the single gg/2 orbit. The effective interaction 
for this orbit has been derived from the CD-Bonn NN 
potential without using any adjustable parameter. This 



approach reproduces very well the excited states of ^^Pd 
observed in a recent experiment, and therefore provides 
confidence in our predictions for ^^Cd and ^^In. 

Aside from the intrinsic interest in employing a realistic 
effective interaction to describe these neutron-deficient 
nuclei, the present work was motivated by the suggestive 
interpretation given in Ref. [l| of the almost equidistant 
levels in the ^^Pd spectrum. It was argued, in fact, that 
this feature may be traced to an isoscalar spin-aligned 
np coupling scheme replacing the normal isovector J = 
pairing which is dominant for like valence particle nuclei. 
We have thus decided to investigate the role played by 
the np interaction in the J = 9, T — channel. 

As regards ^^Cd, we have found that the matrix el- 
ement Vg{np) has a direct influence on energies of the 
high-spin states 12+, 14+ and 16+, making the latter 
isomeric, while for the states with I^ from 0+ to 10+ 
it also substantially affects the structure of the wave 
functions. More precisely, a more attractive matrix el- 
ement leads to a larger content of J = 9 pairs and, 
for states up to /'^ = 8+, to a larger fragmentation 
in terms of the [nn] [pp] basis. It should be men- 
tioned, however, that the large content of J = 9 pairs for 
I^ = 0+, 2+ and 4+ arises also from the significant over- 
lap of the [(np)9{np)9]i component with [{nn)0{pp)I]i 
and [{nn)I{pp)0]i. 

In agreement with previous papers, we have found that 
the J = 9 matrix element plays an important role in 
determining the low-energy spectrum of both ^^Cd and 
^^Pd. However, as shown in some detail for ^^Pd, it does 
not contribute to the energy spacings directly, but rather 
by changing the weights of the isovector contributions 
through the induced fragmentation. So to say, the pair- 
ing force pushed out the door comes back through the 
window. 

In summary, we confirm the relevant role of the 
isoscalar spin-aligned coupling in the low-energy states 
of ^^Pd and ^^Cd. Based on our study and the presently 
available data for ^^Pd, we feel, however, that one can 
hardly speak of a new phase of nuclear matter similar to 
the well-known one induced by the strong pairing corre- 
lations between identical particles. 
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